Generalized Elliott-Yafet theory of electron spin relaxation in metals: the origin of 
the anomalous electron spin life-time in MgB 2 
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The temperature dependence of the electron spin relaxation time in MgE>2 is anomalous as it does 
not follow the temperature dependence of the resistivity above 150 K, it has a maximum around 
400 K, and it decreases for higher temperatures. This violates the well established Elliot- Yafet 
theory of electron spin relaxation in metals. We show that the anomaly occurs when the quasi- 
particle scattering rate (in energy units) becomes comparable to the energy difference between 
the conduction- and a neighboring band. We find that the anomalous behavior is related to the 
unique band structure of MgB2 and the large electron-phonon coupling. The saturating spin-lattice 
relaxation can be regarded as the spin transport analogue of the Ioffe-Regel criterion of electron 
transport. 
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Knowledge of the electron spin-lattice relaxation time, 
Ti, of conduction electrons plays a central role in as- 
sessing the applicability of metals for information pro- 
cessing using electron spins, spintronics T\ is the 
time it takes for the conduction electron spin ensemble 
to relax to its thermal equilibrium magnetization after a 
non-equilibrium magnetization has been induced e.g. by 
conduction electron-spin resonance (CESR) excitation |2j 
or by a spin-polarized current 0. The Elliott-Yafet (EY) 
theory of T% in metals [H, 0| has been well established in 
the past 50 years on various systems such as elemental 
metals @, strongly correlated one-dimensional @, and 
some of the alkali fulleride salt metals. It is based 
on the fact that the spin part of the conduction electron 
wave functions is not a pure Zeeman state but is an ad- 
mixture of the spin up and down states due to spin-orbit 
(SO) coupling. As a result, momentum scattering due to 
phonons or impurities induces electron spin-flip, which 
leads to spin relaxation. Typically every millionth mo- 
mentum scattering is accompanied by the electron spin- 
flip due to the relative weakness of the SO coupling. 
Thus, Ti > t (t being the momentum relaxation time) 
which explains the motivation behind the efforts devoted 
to the spintronics applications of metals. 

A consequence of the EY theory is the so-called Elliott- 
relation, i.e. a proportionality between Ti and r [3j]: 
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Here a is a band structure dependent constant and for 
most elemental metals a w 1..10 (Ref. [5j). L is the 
SO splitting for spin up and down electrons in a valence 



(or unoccupied) band near the conduction band with an 
energy separation of AE. E.g. in sodium, the conduction 
band is 3s derived and the relevant SO state is the 2p 
with AE = 30.6 eV and L = 0.16 eV giving (L/AE) 2 = 
2.7- 10" 5 d. 

The Elliott-relation shows that the temperature depen- 
dent resistivity and CESR line-width are proportional, 
the two being proportional to the inverse of t and Ti, 
respectively. This enabled to test experimentally its va- 
lidity for the above mentioned range of metals. Much as 
the Elliott-relation has been confirmed, it is violated in 
MgB2 as therein the CESR line- width and the resistivity 
are not proportional above 150 K [8(. 

Here, we study this anomaly using MgB 2 samples with 
different B isotopes and impurity concentrations and we 
show that the anomalous effect is indeed intrinsic to 
MgB2. We explain the anomaly with an exact treat- 
ment of the SO scattering of conduction electrons in the 
presence of a nearby band with energy separation AE, 
by extending the Mori-Kawasaki formula developed for 
localized spins to itinerant electrons. The result shows 
that the Elliott-relation breaks down when AE is com- 
parable to h/r. Adrian deduced a similar result with a 
qualitative argument 

The role of AE is disregarded in the EY theory since 
typical values are AE sw 10 eV and H/t — 2nk-BT\ « 
6 meV at T = 100 K and A = 0.1 clcctron-phonon cou- 
pling. We show that the occurrence of the anomaly in 
MgB 2 is related to the unique features in its band struc- 
ture and the large electron-phonon coupling. 

We performed CESR measurements on three kinds of 
fine powder MgB 2 with isotope pure 10 B, 11 B, and nat- 
ural boron (20 % 10 B and 80 % n B). The samples have 
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FIG. 1: Comparison of the temperature dependent CESR 
line-width (■: Mg 11 B2, : Mgl>2 of natural boron) and the 
resistance (solid curve) for Mg lx B2. The two types of data 
overlap in the 40-150 K temperature range. A representative 
error bar is shown. 



slightly different impurity content, shown by the varying 
residual CESR line- width, AB^. The temperature de- 
pendent T\ and the CESR line-width, AB, are related: 
AB = AB + 1/771, where 7 /27r = 28 GHz/T is the elec- 
tron gyromagnetic factor. ESR spectroscopy was done on 
a Bruker X-band spectrometer (center field 0.33 T) in the 
4-700 K temperature range on samples sealed under He 
in quartz tubes. The most important result of the cur- 
rent report, the anomalous temperature dependence of 
AB or T\ , is independent of sample morphology, isotope 
content, or thermal history. AB is also independent of 
the magnetic field, apart from a small change in ABq 
(lo| . Resistance on pellet samples and SQUID magne- 
tometry were studied on the same batch as those used for 
ESR. The RRR > 20 and the sharp (< 0.5 K) supercon- 
ducting transition attest the high quality of the samples. 
Heating the samples in the ESR measurement (about 1 
h duration) to 700 K does not affect the superconducting 
properties as shown by magnetization measurements. 

We reported previously the anomalous temperature de- 
pendence of the CESR line- width in Mg n B2: although 
the line- width follows the resistance for the 40-150 K 
temperature range, it deviates above 150 K and satu- 
rates above 400 K [8| • This was confirmed independently 
our knowledge, this is the only known metal 
where such phenomenon is observed. We extended the 
previous measurement to 700 K and the result is shown in 
Fig. [TJ Interestingly, the CESR line-width does not just 
saturate at high temperatures, as found previously, but 
decreases slightly above 500 K. The result is reversible 
upon cooling with no dependence on the thermal treat- 
ment protocol. The phenomenon is reproduced on several 



samples of different purity and boron isotopes, thus it is 
intrinsic to MgB2. 

We explain the anomalous temperature dependence of 
7i in general before including the specifics of MgB2. The 
Elliott- Yafet theory disregards the magnitude of r and 
takes life-time effects only to lowest order into account 
0, 0]. The extended description involves the Kubo- 
formalism and is based on a two-band model Hamilto- 
nian, H = Hq + Hso, where: 



H = ^ [e»(k) + fvyBs] c^ u s c k ,v, s + i7 Bca tt, 

k,L> ,8 



(2) 



Here — 1 or 2 are the band, s, s' are spin indices, 
L 8iS i is the SO coupling, and B is the magnetic field along 
the z direction. i7 sca tt is responsible for the finite r. The 
SO coupling does not split spin up and down states in 
the same band for a crystal with inversion symmetry, 
however it joins different spin states in the two bands 
[l| . The Hamiltonian in Eq. [2] is essentially the same 
as that considered by Elliott [3j. However, instead of a 
time-dependent perturbation treatment, we calculate T\ 
from the Mori-Kawasaki formula [HI, [lj] : 
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where xo is the static magnetic susceptibility, wl = "fB 
is the Larmor frequency, and Gp P+ (uj) is the Fourier 
transform of 



G 



PP+ (t) = -*e(<)([P(t),P+(0)]) Ho , 



P = [H so ,S + }. 



(4) 



The expectation value in Eq. [4] is evaluated with the 
unperturbed Hamiltonian, Hq. 

Assuming that the two bands are separated by 
AE{k) = ei(fc) — 62 (k) = hAcu(k), a standard calcula- 
tion yields [15j: 
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means Fermi surface averaging, L z {k) 



Li^(k) — Li t i(k), and we neglected wl, which is small 
compared to Auj(k F ). Eq. was previously deduced by 
Adrian using a qualitative argument, which involved an 
effective magnetic field, L/frf, fluctuating with t corre- 
lation time due to the SO coupling 

We approximate Eq. [5] using effective values for the 
band-band energy separation and the SO coupling: 
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FIG. 2: (Color online) Schematics of the spin-lattice relax- 
ation in MgE>2 in the two-band model framework. The ar- 
row thicknesses represent the relaxation rates (not to scale). 
Note that the inter-band momentum scattering rate is larger 
than the spin-lattice relaxation rates, therefore there is a spin 
transfer between the two types of bands. 
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This result returns the Elliott-relation when rAw e ff ^ 1 
and gives a decreasing spin relaxation rate with increas- 
ing r _1 when rAu; c ff < 1, thus it can be regarded as 
a generalization of the Elliott- Yafet theory. In the fol- 
lowing, we show that it describes the spin relaxation in 
MgB 2 . 

Electronic properties of MgB2 are described by the 
so-called two-band model meaning that the conduction 
bands related to the boron a and 7r bonds have different 
electron-phonon couplings, different affinity to defects, 
and that the inter-band momentum scattering is weaker 
than the intra-band ones As a result, the conduc- 

tivity is given by a parallel resistor formula [l6j], i-e. the 
band with longer r dominates the transport. In contrast, 
the CESR spin relaxation is dominated by the band with 
shorter T\ . Although the inter-band momentum scatter- 
ing time, is longer than the intra-band momentum 
scattering times, r CT and r w , it is still much shorter than 
T\ . Thus an electron with a given spin state is scattered 
back and forth between the two types of bands several 
times before flipping its spin, which is depicted in Fig. [2] 
The overall 1/T\ is the average of the spin-lattice relax- 
ation rates weighted by the relative DOS on the a and 7r 
bands, N„ = 0.56 and N a = 0.44 0: 
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In Fig 
Refs. m 



3, we show the band structure of MgB2 from 
19j near the Fermi energy. Two boron a and 




FIG. 3: (Color online). Band structure of MgB2 near the 
Fermi energy after Refs. [IH, [ljj. Two of the a bands (black) 
cross the Fermi surface near each other in the vicinity of the 
r and A points, whereas ir bands (red) are separated from 
other bands with a larger optical gap at the crossing. We also 
show the dispersion with 8 times larger wave- vector resolution 
around the T points with vertical arrows for possible AE a 
values. 



two 7r bands cross the Fermi energy such that the it bands 
are well separated from other bands with AE^ > 2 eV 
whereas the two a bands are close to each other and 
AE a « 0.2 eV. Based on the above theory and Eq. [BJ, 
we conclude that I\ follows the EY mechanism for the 
7r bands, whereas it is described the by the novel mech- 
anism for the a bands. With this in mind and the two 
band model result of Eq. [3 we describe the CESR line- 
width with: 



AB = AB 
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where we introduced the band index for the parameters. 
The momentum relaxation times are calculated using the 
Debye-model and assuming clean samples, i.e. zero resid- 
ual scattering: 
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where n = a, ir, tojy is the Debye frequency, and Xtr,n are 
the transport electron-phonon couplings from Ref. [l6j |. 
which contain both intra- and inter-band scattering. 
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TABLE I: Parameters used to calculate the CESR line-width 
in MgB2. The given standard deviations indicate the free 
parameters of the fit. 



200 



400 
T(K) 



600 



FIG. 4: Measured (symbols) and calculated (solid curves) 
CESR line-width in MgB 2 with n B and 10 B. Note the larger 
residual line-width in the latter sample. Dashed curves show 
the contributions to the line-width from the a and n bands 
for the 11 B. Dotted curve shows a calculation for the the 11 B 
sample assuming 1/Ti is due to a bands only. 



In Fig. H, wc show the CESR line-width for Mg n B 2 
and Mg 10 B 2 between 40 and 700 K and the calculated 
line- width using Eq. [H with parameters in Table [U ob- 
tained from a fit. Results on the natural boron sample 
are identical to the data on the Mg n B2 within exper- 
imental error and are not shown. The larger residual 
line-width in the 10 B (AB Q = 2 mT) than in the n B 
sample (AB = 1 mT) is related to a larger defect con- 
centration in the starting boron, the preparation method 
and the starting Mg being identical. Apart from this, 
the only difference between the two samples are the dif- 
ferent Debye temperature, Od- The calculated CESR 
line-width (solid curves) reproduces well the experimen- 
tal data with the parameters in Table U The dotted 
curve in Fig. 2] is a calculation assuming that relaxation 
is given by the a bands alone, which accounts relatively 
well for the data with three free parameters (L CT , AE e g tCr , 
and AB Q ). However, it fails to reproduce the slope of AB 
at higher temperatures, which shows the need to include 
relaxation due to the n bands. 

The determination of AE C B,a ~ 0.2 eV is robust as it 
is given by the temperature where the maximal AB is 
attained and its value is close to values expected from 
the band structure (arrows in Fig. [3]). Knowledge of 
AE D f[^ allows to determine the SO splitting indepen- 
dently, -Lcff.cr = 0.64 meV, as usually only the L/AE 
ratio is known. The SO splitting for the atomic boron 
2p orbital is L = 0.23 meV (Ref. 0), which is in a rea- 
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sonable agreement with the experimental value. AE^ 
was fixed to 2 eV which affects T G ff j7r as these are not 
independent. 

The isotope effect on O p is 10 e D / 11 OD = 1.04, that 
is close to the expected a/H/10 ratio. The Od values 
are in agreement with the 440.. 1050 K values in the lit- 
erature, which scatter depending on the experimental 



method 2fj, 21 1. We note that the model could be im- 
proved by including the Einstein model of phonons or by 
an exact treatment of the band structure dependent SO 
coupling and band-band separation. 

Finally, we note that the maximum of 1/Ti occurs 
when tAuj rj 1. This coincides with the Ioffe-Regel crite- 
rion for the electron transport [23l | when the band-band 
separation is comparable to the bandwidth, w, e.g. in 
narrow band metals. For MgB2, w m 10 eV [18J therefore 
saturation of the CESR line- width is not accompanied by 
a saturation of electrical resistivity. 

In conclusion, we explained the anomalous spin-lattice 
relaxation in MgB 2 by extending the Elliott- Yafet theory 
to the case of rapid momentum scattering and near lying 
bands. The anomaly does not occur in conventional met- 
als, which have small electron-phonon coupling and well 
separated bands. A similar phenomenon, the so-called 
Dyakonov-Perel relaxation [l| , occurs for semiconductors 
without inversion symmetry, although its physical origin 
is different. The band structure of some of the other di- 
borides in e.g. BeB2 and CaB 2 predicts [191 ] similar phe- 
nomena but conventional spin relaxation in A1B 2 , ScB 2 , 
and YB 2 . We also predict that the described effect is 
sensitive to pressure since this shifts the a bands [24j . 
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